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ABSTRACT 

The  purpose  of  this  paper  is  to  describe  and  illustrate  a  program, 
written  for  the  IBM  7090,  which  defines  a  model  for  a  finite  collection 
of  algebras  for  the  computer.  It  is  shown  that  the  program  contains  a 
structure  broad  enough  in  scope  to  allow  one  to  perform  operations  on 
such  diverse  mathematical  concepts  as  differential  equations,  infinite 
series,  and  differential  forms  in  a  simple  yet  comprehensive  manner, 
while  also  serving  as  a  foundation  upon  which  a  variety  of  higher- 
level  symbolic  manipulation  languages  can  be  developed. 
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Introduction 

It  has  become  increasingly  evident  in  recent  years  that  the  intri¬ 
cate  numerical  procedures  that  have  been  designed  for  the  computer  have 
not  been  adequate  for  handling  mathematical  and  physical  problems  requir¬ 
ing  extensive  symbolic  analysis;  what  has  been  needed  is  the  development 
of  highly  sophisticated  symbolic  techniques  treating  various  analytical 
problems  in  a  variety  of  ways.  Studies  have  been  made  regarding  this 
problem,  but  because  of  the  many  inherent  difficulties  involved,  much' 
confusion  has  arisen  regarding  not  only  the  nature  of  the  computational 
procedures  that  could  best  be  handled  using  the  computer,  but  the  structure 
that  would  be  necessary  for  the  development  of  any  higher-level  language 
which  would  be  designed  to  manipulate  such  data.  In  the  present  paper, 
we  shall  attempt  to  rectify  this  situation  in  part  by  defining  a  model 
for  a  finite  collection  of  (associative)  algebras.  It  will  then  be 
seen  that  the  program,  called  ALGEBRA,  will  provide  us  with  a  structure 
broad  enough  in  scope  to  allow  us  to  perform  operations  on  such  diverse 
mathematical  concepts  as  differential  equations,  infinite  series,  and 
differential  forms  in  a  fairly  simple  yet  comprehensive  manner.  More¬ 
over,  the  model  will  serve  as  a  foundation  upon  which  a  variety  of 
higher-level  languages  dealing  with  symbolic  manipulation  can  be  con¬ 
structed,  depending  only  on  the  types  of  mathematical  computational 
procedures  desired. 

Before  we  proceed,  however,  several  remarks  should  be  made  regard¬ 
ing  the  type  of  programming  language  that  is  appropriate  for  developing 


1 


symbolic  computational  procedures  of  any  type.  In  general,  assembly 
languages  are  unsuitable  because  of  the  complexity  of  even  the  most 
elementary  operations  on  complex  data  structures  that  are  generally 
required,  whereas  the  algebraic  languages,  such  as  FORTRAN  and  ALGOL, 
are  excluded  because of  their  relatively  primitive  storage  allocation 
capabilities.  The  list  processing  languages,  however,  are  quite  adequate 
for  this  purpose  because  of  their  ability  to  store  and  manipulate  exceed¬ 
ingly  complex  data  structures.  Of  the  systems  in  this  category,  the 
LISP  1.5  language  (see  reference  1  or  6)  was  used  in  the  development  of 
ALGEBRA  because  of  its  simplicity  and  generality. 

Initial  Formulation  of  the  Model 

Let  ft  =  {R1,...,Ra}  be  a  set  of  commutative  rings  and  h  =  [Ax  , . 
a  collection  of  algebras,  each  of  which  is  operated  on  by  a  ring  in  R 
(see  reference  3).  It  will  be  assumed  that  the  rings  and  algebras  have 
a  multiplicative  identity  "1"  and  that  each  algebra  K,  is  an  algebra  of 
some  module  generated  by  a  recursive  set  -1  For  convenience  of  refer¬ 
ence,  the  algebras  will  be  indexed  by  the  variable  v  and  the  rings  by 
the  variable  r.  Moreover,  for  each  ring  R.  of  R,  a  computable  function 
s  -•  scalar,,  (s)  will  be  assumed  to  exist  which  evaluates  expressions  s 
of  the  forms  aa  +  . . .  +  Sp ,  ax  *  •  •  aj, ,  and  aq  ,  where  a,a1>...,a?  are 
scalars  (i.e.,  elements  of  the  ring)  and.  q  is  a  non-negative  integer. 

In  general,  the  definition  of  this  mapping  will  be  extended  to  evaluate 

1A  set  A  will  be  said  to  be  recursive  if  for  any  expression  x,  the 
predicate  x  e A  is  computable  (see  reference  4  or  5). 
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re  elaborate  expressions  containing  the  scalars,  but  this  will  neces- 
rily  vary  with  the  structure  of  R„  and  how  we  wish  to  use  it.  In  any 
se,  the  function  will  always  serve  as  the  sole  criterion  for  whatever 
sumptions  are  being  made  regarding  this  ring. 

Consider  now  an  algebra  AT  defined  over  the  ring  R.  ,  the  elements 
which  will  be  called  vectors.  We  shall  find  it  useful  for  a  variety 
reasons,  including  input-output  and  uniqueness  considerations,  to 
sume  that  there  exists  a  linear  ordering  relation  of  the  set  . 
may  also  wish  to  presume  that  the  algebra  is  either  free  or  semi- 
mple,  though  this  will  not  be  necessary.  It  will  be  convenient,  how- 
er,  to  be  able  to  represent  the  vectors  in  various  forms,  but  this 
so  occurs  whenever  the  theory  is  developed  from  a  mathematical  view- 
•  int. 


Now  any  non-zero  vector  x  in  A.  can  always  be  written  -in  the  form 
*  ca  (tpi )  ^  1  *  •  •  (cpi  r  l  +  ...  +  cn  (cpy)AnJ  '  *  *  (epp  n  )  "  r  n  ,  where  each 
/  0  is  a  scalar,  \t,  is  a  positive  integer,  and  cp!  eX,.  W  {l } .  If  the 


gebra  A.  is  commutative,  it1  (will  be  presumed  that  the  terms  of  the 

A.  ■  ?  !  ■ 

•oduct  (tpj)  1 1  . .  •  (cp.1,  )  '• r  i  are  ordered  so  that  cpj  <  cp^  <  ...  <  cpi  ; 
'.herwise,- their  ordering  wiljL  have  to  be  specified  in  the  hypothesis 
.ven  concerning  the  algebra.  In  any  case,  with  the  further  assumption 
rnt  (cpl)*11  ••'(cpx)Xl1  <  •••  <  (cpr1)*’lr*  •••  (cpn  )  nra>  it  follows  that 

n 

ssociated  with  every  vector  is  an  expression  of  the  form  x  =  ((ci,(cpi ) , 


»  CtPr  x  )  1  r  1)  , . . . ,  (cn/r{)  nl  ,  . . . ,  (cp*r  )  "‘a)).  This  form  will  be  called 
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the  canonical  representation  of  x  and  will  be  denoted  by  partfx.X.,). 

For  completeness,  we  shall  also  require  partfOjX.,)  to  have  the  value 

«0,1». 

Whereas  the  mathematical  form  of  a  vector  may  or  may  not  be  useful 
for  an  analyst,  depending  on  the  nature  of  the  problem  under  considera¬ 
tion,  the  form  will  always  be  found  to  be  awkward  for  a  computer  be¬ 
cause  of  the  partitioning  that  necessarily  will  ensue.  Consequently, 
the  usual  procedure  for  manipulating  these  elements  is  to  first 
partition  them  into  their  canonical  representations,  use  these  repre¬ 
sentations  for  carrying  out  the  operations  that  are  needed,  and  then 
return  the  answers  in  their  mathematical  forms.  In  so  far  as  the 
uniqueness  of  the  canonical  forms  is  concerned,  and  this  will  be 
discussed  in  detail  later,  let  it  suffice  to  say  for  the  moment  that 

this  will  most  certainly  depend  on,  among  other  things,  the  assumptions 

. .  .*  * 

made  regarding  tne  basic  structure  of  A, . 

We  now  consider  the  addition  and  multiplication  operations  for  the 
vectors  of  A.,,  which,  for  convenience,  will  be  denoted  by  the  symbols  ® 
and  *  respectively.  The  addition  function  is  particularly  easy  to  de¬ 
scribe,  being  defined  as  follows: 

1 .  For  any  vector  x,  x#0  =  0$x=x 

2.  If  x  is  a  non- zero  vector  having  the  canonical  form 

<(a,<1,...,  y  is  an  arbitrary  non-zero  vector 

having  the  canonical  form  (  (bx  ,  (cpj  )  ^  1 ,  .  . . ,  (cpj  )  r  1 )  ,  .  .  . , 


,  ^'n  r  n  \  \  e  _  /  -i  \  A.»  1  ,  1  v  ^1 


(b»»(<PS>  »•••*(&')'*'*)),  and  ?J  =  (cp‘)Alx _ <cp^  t)  lr‘ 

(i  ~  1  j  •  jii)  y  then 


+  •••  +  h1§!  +  aY,1  ...  Y^  +  bt  +  \5i  +  i  +  •••  +  bn§n 

if  5,  <  rf1  ...  y“=  <  s„. 


J  b^l  +  *’•  +  bi-i5l-i  +  scalar,  (a+b,)*^  +  b1+1§1  +  l  +  bn  §n 
if  =  Yj1  ...  Y3  =  and  scalar,  (a+bj)  ^  0 

bi§i  +  •••  +  +  bi  +  i§1  +  i  +  ...  +  bn§r 

if  §1  =  ...  and  scalar,,  (a+bj)  =  0 

3.  If  x  is  a  vector  having  the  canonical  form  ((ax ,  (Y;)^1  , . . . , 
j)  ),...,  (aTl  (Yi )  ,  •  •  •  >  (^j  _  )^a  *  * )) ,  y  is  an  arbitrary  vector, 

(y{  )^!  1  . . .  (Ys  j  )^!s  1 ,  then  x  ®  y  ^  (a2§2$(...»  (a3  £,fcy)  . . .))  . 

illarly,  using  the  associative  and  distributive  properties  of  the 
;ration  *,  it  is  clear  that  the  definition  of  x  *  y  reduces  to  the  case 
jre  x  has  the  canonical  form  ( (a, Y^1 , . . . ,Y,a ))  and  y  has  the  canonical 
rm  ( (b,c^L  1  ,  ...,cftN)).  However,  since  the  definition  of  this  product 
:ies  with  the  algebra  under  consideration,  it  is  evident  that  it  will 
/e  to  be  given  by  a  function  x,y  -  vectorv (x,y).  This  map  will  always 
assumed  to  be  computable,  of  course.  . 

1MPLE: 

Let  R  be  the  standard  fixed  point  arithmetic  indexed  by  the  value 
»  0,  so  that  the  function  s  ->  scalar0(s)  has  as  its  value  the  sum 
roduct)  of  the  numbers  a1,...,ap  whenever  s  «*  ax  +  ...  +  a„ 

“  ai'^Ap)  and  the  exponentiation  of  the  number  a  to  the  non-negative 


{»;*  '  -  <**•  " 


•  <  *  1  "**  <■*»«»*■  ' 


integral  power  n  whenever  s  =  ar‘  .  If  A  is  the.  polynomial  algebra  of 
the  set  X  =  {x)  over  the  ring  R  (see  reference  2),  it  is  then  evident 
that  the  mapping  X,p  -♦  vector  (X,p)  can  be  easily  defined  by  assigning 
the  value  scalar0(a*b)xSCa^arc>  ^+-0  sealary  (k+£)  f  0  and  scalar^ (a-b) 
if  scalar0  (k+£)  -  0,  where  X  =  axk  and  p.  =  bx  .  Consequently,  it 
follows  that  the  expressions  K  +  L  and  K*L  easily  reduce  to  the  respec¬ 
tive  values  7  +  Ax  +  3xs  and  12  +  16x  +  9x2  +  12x3  for  K  =  3  +  4x  and 
L  =  4  +  3X2  . 

Functions 

Let  K,  be  an  algebra  in  h  defined  over  a  ring  Rr .  A  collection  of 
vector-valued  functions  (with  vectors  and/or  scalars  as  arguments)  is 
generally  associated  with  this  algebra,  some  resulting  directly 
from  the  algebraic  structure  of  Av  and  others  depending  on  whatever 
further  restrictions  are  placed  on.A, ..  In  any  case,  since  the  choice 
of  these  maps  will  necessarily  vary  with  the  algebra  under  consideratic 
and  how  it  is  to  be  used,  it  follows  that  the  model  ALGEBRA  will  have 
to  contain  a  general  procedure  for  storing  and  manipulating  any 
assumptions  that  will  be  required. 

Before  the  vector-valued  functions  are  considered,  however,  certair 
remarks  should  be  made  regarding  the  structure  of  the  algebra.  In  par¬ 
ticular,  although  it  is  true  that  in  certain  cases  we  shall  be  interest 
in  Ay  for  its  algebraic  properties  only,  it  is  also  clear  that  at  other 
times  the  topological  aspects  of  the  algebra  will  be  of  dominant  im¬ 
portance.  For  example,  the  assumption  that  the  algebra  is  a  Banach  spe 
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is  a  fairly  common  requirement,  in  which  case,  we  would  most  certainly 
be  interested  in  the  bounded  linear  operators  of  this  space.  In  fact, 
we  might  even  be  interested  in  the  vector-valued  Borel  functions.  In 
many  cases,  however,  the  functions  that  we  shall  be  interested  in, 
whether  they  are  continuous  or  not,  will  not  be  computable.  And  for  that 
matter,  we  might  be  interested  in  using  a  computable  function  for  its 
representation  only,  and  not  its  value. 

Heretofore,  it  has  been  assumed  that  A.  is  an  algebra  of  some  module 
K,  ,  where  M^.  is  generated  by  the  set  X..  .  If,  however,  we  wish  for  our 
algebra  to  have  a  topological  structure,  this  assumption  may  be  relaxed 
by  requiring  that  X.  be  a  set  that  is  dense  in  ML .  The  set  X,  can  be 
modified  in  other  ways  also.  In  particular,  consider  a  function  / 
associated  with  A,.  Then  for  any  expression1  s  of  the  form  /(x.  ,  .  .  .  , x.  )  > 
where  each  x<  is  a  vector  or  scalar,  if  it  isn't  desirable  to  compute 
s  for  the  arguments  under  consideration,  or  it  is  not  possible  to  do  so, 
the  expression  will  be  regarded  as  being  irreducible  relative  to  the 
hypotheses  being  made,  and  hence  will  be  treated  as  an  element  of  X^.  . 


1  An  expression  in  \  is  defined  to  be  (1)  a  scalar  or  vector  in  A.,  ,  (2) 
a  representation  of  one  of  the  forms  xz  +  ...  x.  or  X)  . . .  x.  ,  where 
each  x.  is  an  expression,  (3)  a  representation  of  the  form  x-  ,.  where 
x  is  an  expression  and  n  is  a  positive  integer,  or  (4)  a  representation  of 
the  form  /(xx , . . . , xs ) ,  where  /  is  an  associated  mapping  and  each  x{  is  an 
expression. 
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In  general,  if  no  assumptions  are  made  concerning  the  function,  we 

have  no  choice  but  to  regard  it  as  being  irreducible  whenever  it  is 

used.  If,  on  the  other  hand,  /  is  to  satisfy  certain  hypotheses,  a 

function  h  is  associated  with  /.  This  mapping,  called  the  axiom  or 

hypothesis  function  of  /,  acts  as  a  partial  evaluation  function  for  /, 

stating  which  properties  of  the  function  are  to  be  assumed,  and  how 

they  will  be  used.  In  particular,  if  /  has  the  property  e. (xj , . . . , x. ) 

whenever  the  condition  pj  (xj , . . . ,xB )  is  fulfilled  (for  i.=  l,2,...,n), 

* 

then  h  is  the  function  of  the  variables  x1,...,x-,  having  for  its  value 
et (xj  »...»Xb)  if  the  predicate  p.  (xx , . . . ,xa)  is  satisfied,  and  the 
irreducible  expression  /(x1,...,x=)  otherwise. 

Let  /i >•••»/„  be  a  collection  of  functions  which  are  to  be  associat 
with  A,  and  which  will  satisfy  the  axiom  functions  hj,...^,.,  respectiv 
These  functions  and  their  properties  are  referred  to  in  ALGEBRA  by  a 
variable,  called  AXIOM,  whose  value  is  the  list  (/i»h1,/2,h2,i..,/n,h!;) 
If  no  assumptions  are  to  be  made,  AXIOM  will  have  as  its  value  the 
empty  set  0.  In  any  case,  whenever  a  vector-value  expression  / (x1  , . . . , 

I 

is  encountered,  iJ:  it  is  not  a  sum,  product,  or  exponential  expression 
of  the  form  xy  where  x  is  a  vector  and  y  is  a  non-negative  integer, 
ALGEBRA  first  checks  to  see  if  /  is  a  member  of  the  set  AXIOM.  If  it 
is  found  to  be  in  AXIOM,  the  expression  is  evaluated  using  the  corre¬ 
sponding  axiom  mapping  of  /;  otherwise,  s  is  regarded  as  being  irreduci 
and  hence,  a  member  of  X,  . 
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EXAMPLE: 


Let  A  be  the  polynomial  algebra  defined  in  the  example  on  Page  5 
and  assume  that  the  trigonometric  functions  sin  and  cos  are  associated 
with  A.  If  they  are  assigned  the  respective  axiom  functions  hsin  and 
hcos,  where 


hsin(s)  =  < 


"0 _ : _ _ _  _ if  s  =  .0  .  . _ _ _ . . . 

(-l)*hsin[ (-l)*s]  if  s  =  n  or  n§  (n=-l,-2,...  and  §  a  vector) 

hsin(§)*hcos[ (n-l)§]©hcos (§)*hsin[ (n-l)§] 

if  s  =  n§  (n=2,3,...  and  §  a  vector) 
hsin(§x  )*hcos  (?2+. .  .+§„  )©hcos  )*hsin(§2  +  *  •  •+§n  ) 

if  s  =  §!+.■.+§„  (each  a  vector) 
sin(s)  otherwise 


and 

”l  if  s  =  0 

hcos[(-l)*s]  if  s  =  n  or  n|  (n=-l , -2 ,  . . .  and  §  a  vector) 

hcos  (0*hcos[(n- 1  )§]©(-!  )*hs  in  (§)*hsin[(n-l)§}  . . . 

if  s  =  n§  (n=2,3,...  and  §  a  vector) 
hcos (§! )*hcos (§2+. . .+En)©(-l)*hsin(§1)*hsin(|2+. . .+§„) 

if  s  =  §i+...+§n  (each  5,  a  vector) 

^cos(s)  otherwise 

T  _  s  i 

it  then  follows  that  any  expression  s  of  the  form  ^ at [sin(p1 x) ] 


hcos(s)  =  < 


[cos (qt x) ]" 1  (aj,Pj,  and  q{  are  arbitrary  integers  and  mt ,  n;  =  0,1,2,...) 
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easily  reduces  to  a  vector  of  the  form  £  bj[sin(x)]  ,i[cos(x)] 

Moreover,  if  it  is  also  assumed  that  sin2 (§)  +  cos2(§)  =  1  for  every 
vector  §  in  A  (i.e.,  the  mapping  \,(ji  -*  vector  (A, p.)  is  modified  by  letting 
vectorx [cos (§) ,cos(§)]  =  1  -  [sin(§)]2>,  then  s  completely  collapses  to 
an  expression  of  the  form  £  o-r  [sin(x) ^  +  £  3,  [sinCx)]^5  cos<x)  .  Thus, 
we  see  that  the  trigonometric  polynomials  over  the  scalar  ring  R  are 
easily  manipulated  whenever  this  algebra  is  used. 

Definition  of  the  ALGEBRA  System 

In  the  previous  two  sections  a  procedure  was  developed  for  first 
formulating  the  basic  structure  of  an  algebra  Ay,  and  then  associating 

with  this  algebra  an  arbitrary  collection  of  vector-valued  maps,  along 
with  any  associated  assumptions  that  were  required.  The  generating  set 

X,  of  this  algebra  was  seen  to  be  well-defined,  consisting  entirely  of 
explicitely  defined  vectors  and  irreducible  expressions  of  the  form 
/(xx , . . . ,x„) .  Since  these  elements  were  evaluated  to  themselves,  it 
followed  recursively  (by  the  definitions  of  the  ®,*,  and  axiom  mappings) 
that  every  combination  s  of  expressions  of  these  elements  would  reduce  to 
an  unique  value,  the  form  of  which  would  vary  only  when  the  underlying  struc¬ 
ture  was  changed.  We  now  define  the  ALGEBRA  system  to'  be  this  reduction 
procedure,  in  which  case  it  is  then  representable  as  a  function  of 

the  form  s,  X,  ,r,v, AXIOM  -  ALGE BRA ( s, X.,  ,r,v, AXIOM)  .  - — - — - — 

Let  s  be  an  expression  in  A,  having  the  value  ALGEBRACs,^  ,r,v, 

AXIOM)  =  x.  Since  x  and  s  obviously  represent  the  same  element  of  the 
algebra,  it  is  clear  that  the  mapping  defines  an  equivalence  relation 
on  the  set  of  expressions  of  A,,  requiring  any  two,  say  sa  and  s2 ,  to 
be  equivalent  whenever  ALGEBRA^  ,  X^  ,r,v, AXIOM)  =  ALGEBRA(sa  jX^  ,r ,v, 

AXIOM).  Consequently,  since  the  hypotheses  concerning  the  structure 
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of  the  algebra  have  been  imbedded  in  the  functions  a  -*  scalar,,  (a)  and 
A,|j,  -*  vectorv  (\,p.)  and  in  the  variable  AXIOM,  it  follows  that  the 
function  ALGEBRA  does  indeed  characterize  the  structure  of  A,  as  we 
have  defined  it.  Moreover,  since  the  values  of  this  mapping  are  unique, 
they  serve,  in  effect,  as  representations  for  their  equivalence  classes. 
Thus,  if  the  definition  of  the  canonical  form  of  a  vector  s  is  modified 
to  be  that  of  the  vector  ALGEBRA ( s, Xy ,r,v, AXIOM) ,  re  then  note 
that  the  canonical  form  of  a  vector  will  always  be  unique. 

In  reflecting  on  the  structure  of  the  mapping  ALGEBRA,  it  should 
first  be  emphasized  that  all  we  have  succeeded  in  doing  is  implementing 
a  model  for  a  finite  collection  of  algebras  on  the  computer;  no 
techniques  have  been  developed  for  proving  theorems  about  these  algebras. 
The  question  of  theorem-proving  in  general  is  a  distinct  problem,  re¬ 
quiring  the  formulation  of  various  techniques  for  handling  the  learning 
and  choosing  procedures  that  are  necessarily  involved.  Secondly,  many 
of  the  restrictions  that  were  placed  on  the  collection  of  sets  h  were 
not  really  necessary.  It  is  evident,  for  example,  that  ALGEBRA  could 
have  easily  been  defined  to  handle  a  collection  of  groups  or  monoids 
instead.  In  any  case,  having  actually  formed  a  simple  but  comprehensive 
model  of  a  collection  of  mathematical  systems,  each  of  which  contains  a 
fairly  complex  underlying  structure,  we  have,  in  effect,  provided  our¬ 
selves  with  an  extremely  powerful  tool  for  treating  a  vast  array  of 
problems  in  a  systematic  manner  which  could  not  even  have  been  considered 
otherwise . 
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v’lr.sj-wo 


.wo,g^acnjic ..'A  j 


Let  R  be  an  arbitrary  ring  indexed  by  the  value  v  =  v0  and  described 
by  the  mapping  s  -*  scalarr  (s),  X  a  recursive  collection  of  non-numerical 

O 

atoms  or  indexed  variables  designated  by  the  value  v  =  v0  (see  reference  1) 
and  A  an  extension  of  the  polynomial  algebra  defined  on  the  set  of 
variables  X  over  the  ring  R  (see  reference  2),  allowing  in  addition  any 


appropriate  associated  functions  that  are  desired.  It  then  follows  that 

A  will  be  well-defined  whenever  the  product  mapping  x,y  -  vectorv  (x,y) 

0 

is  evaluated.  But  this  function  is  easily  defined,  having  for  the  non¬ 
zero  vectors  x  =  acpp  . . .  .<pa*  and  y  =  bYxx  ...YBn  the  value 
scalar,  (a*b)cpP*(cp^2*(. .  .*(cp*  "*(Y* 1  ••  *Y*n  ))...)) , 


where 


cpr  ■*  (Y^ 1  . .  .Y^n)  » 


Y*1...^1  cprY^...Y*n  if  <  <p  <  Ya 

...Y*n  if  cp  =  Ya  and 
ALGEBRA  (r+sa ,  x ,  r0  ,  v0  ,  AXIOM)  =  0 
scalarrQ  (cpr  +  3Qr)'i^1  .••Y^i1Y(^+1  ...Yn°  if 
cp  =  Ya,  cp  and  r  +  Sq,  are  scalars. 


and  cp 


r+»  a 


is  an  expression  in  our 


ring  of  scalars 

YS11...^I1cpr+SQrYsQ,^1...Y,nn"  if  cp  =  Ya  and 
cpr  +  SQf,  is  irreducible. 

ALGEBRA(cpr  +  *a,X,r0,v0,  AXIOM)* 

(Vs  1  /O'-!  ysQ'  +  1  v'M  if  CD  =  Y 

ui  Ya  +  i  •  •  •  ”n  '  t  Ja 

and  cDr  +  #0f  is  reducible 


We  now  consider  several  useful  applications  of  this  algebra* 


1)  The  SIMPLIFY  Algebra 

This  algebra  is  the  algebra  most  generally  used  for  manipulating 
symbolic  mathematical  expressions,  having  for  its  scalar  ring  the  stand¬ 
ard  fixed  and  floating  point  arithmetic  and  for  its  generator  set  X  the 
collection  of  non-numerical  atoms  which  are  not  being  used  as  vectors 
for  other  algebras.  For  convenience,  we  shall  denote  the  value  of  its 

implementation  mapping  s  -*  ALGEBRA(s,X,r0 ,v0 , AXIOM)  by  SIMFLIFY(s). 

ds 

Assume  that  the  standard  differentiation  function  diff:s,x  -* 
is  now  associated  with  SIMPLIFY  (see  reference  2).  If  AXIOM  *  0  (the 
empty  set),  it  is  then  clear  that  although  SIMPLIFY (K  +  x)  -  5x  +  5y  +  3xy 
for  the  vector  K  =  4x  +  5y  +  3xy  (where  x  <  y  <  xy),  since  no  assumptions 
are  being  made  regarding  the  differentiation  mapping,  SIMPLIFY  [diff(K,x)] 
has  for  its  value  diff(4x  +  5y  +  3xy,x),  considering  this  expression  as 
being  irreducible  and  hence  an  element  of  X.  If,  on  the  other  hand,  h(s,x) 
is  the  standard  axiom  mapping  for  diff1  and  AXIOM  =  (diff,k),  then 
SIMPLIFY  [diff(K,x)]  would  have  as  its  value  h(K,x)  =  4  +  3y.  Thus,  if 
the  two  trigonometric  functions  sin  and  cos  are  also  associated  with 
SIMPLIFY  and  h[sin(x),x]  is  defined  to  have  the  value  cos(x),  it  is  evident 
that  whereas  this  algebra  is  appropriate  for  computing  expressions  of  the 
form  cos(x)  +  diff [sin(x) ,x] ,  it  is  totally  inadequate  for  handling  either 
differential  equations  or  polynomials  with  symbolic  coefficients. 


1  h(s,x)  is  defined  for  any  expression  s  and  any  atomic  vector  x  by  requiring 
that  (1)  h(s,x)  =  0  whenever  s  is  a  scalar  or  atomic  vector  different  from  x, 
(2)  h(x,x)  =  1,  (3)  h(s,x)  =  h(si,x)  +  ...  +  h(sp,x)  for  s  =  sx  +  ...  +  sp, 
and  (4)  h(s,x)  =  h(s1  ,x)ss  *  *  *sp  +  ...  +  h(sp,x)  for  s  =  Sjf-s,. 
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2)  The  MULT  and  MULTILINEAR  Algebras 

The  MULT  algebra  uses  as  its  scalar  ring  the  SIMPLIFY  algebra  and 

as  its  generator  set  a  finite  collection  {x1}...Xb]  of  non-numerical 

atoms  or  indexed  variables.  The  standard  differentiation  mapping 
ds 

diff:s,x  -*  is  again  associated  with  this  algebra,  being  defined 
as  before  whenever  x  =  Xj  for  some  xt .  The  definition  of  diff,  however, 
can  be  extended  to  also  cover  the  case  where  x  is  a  non-numerical  atom 
in  the  ring  SIMPLIFY.  This  is  done  by  requiring  that:  (1)  if  s  is  a 
scalar,  then  diff(s,x)  will  be  computed  as  an  expression  in  SIMPLIFY,  (2) 
diff(s,x)  will  be  regarded  as  being  irreducible  for  any  irreducible 
vector  s  in  MULT,  so  that  in  effect,  each  xt  is  a  function  of  x,  (3) 
diff(s,x)  =  diff(s1,x)  +  ...  +  diff(sp,x)  for  any  s  =  Si  +  ...  +  sp , 
and  (4)  diff(s,x)  =  dif f (sx ,x) s2 • • *sp  +  ...  +  Sj  •••  sp_j  diff(sp,x) 
for  any  3=8!  • • •  sp .  Noting  then  that  the  expression  xxt  + 
diff [sin (x)xt ,x]  reduces  to  [x  +  cos(x)]x1  +  sin(x)diff (xt ,x) ,  it  is 
evident  that  any  differential  equation  with  numerical  coefficients  is 
easily  handled  by  this  algebra. 

In  order  to  manipulate  differential  equations  with  symbolic 
coefficients,  the  MULTILINEAR  algebra  must  be  usdd.  This  algebra  is 
a  modification  of  the  MULT  algebra,  having  as  its  scalar  ring  the 
MULT  algebra  (with  generating  set  {x1,..,xa})  and  as  its  generator  set 
a  finite  sequence  [yi,...,y„}  of  non-numerical  atoms  or  indexed  variables 

each  of  which  is  assumed  to  be  a  function  of  the  variables  [x1,...,xB}. 
The  definition  of  the  differentiation  mapping  diff(s,x)  is  tl  an  modeled  a: 

that  of  the  MULT  algebra,  allowing  x  to  be  either  y.  or  yi  and  requiring 
that  the  expression  diff^Xj)  be  irreducible  for  any  irreducible 


14 


expression  s  .in  the  algebra.  Thus,  it  would  folljow  that  the  expression 
diff  [axj  sin(y_,  )  ,xt  ]  +  di£f[yk,yj]  would  reduce  to  a*sin(yj)  + 
axtcos(yj)diff (y^ ,x()  for  j  f  k. 

Series  Algebras 

Let  R  be  a  commutative  ring  with  multiplicative  identity  1,  Z  the 
ring  of  integers,  and  A  the  collection  of  all  functions  f:Z  —  R  such 
that  for  some  integer  n, ,  f(n)  =  0  for  all  n  <  n,  (the  value  of  nf  will 
necessarily  vary  with  the  choice  of  f).  Denoting  the  value  of  f(n)  also 
by  f n ,  it  then  is  clear  that  A  is  an  algebra  over  R  whenever  (af)„  =  afr , 

(f  +  g\  =  f.  +  g.  ,  Md  (fg)r  =  E  fjg,  for  any  scalar  a  and  any  vectors 

4  “  '  i  +  J  “T1 

f  and  g.  This  algebra  is,  of  course,  isomorphic  to  the  algebra  of  series 

CD 

of  the  form  E  anx"  for  some  variable  x. 
r.fck 

Regarding  n  now  as  a  fixed  variable,  represent  the  elements  of  A  by 

i  •  •  ; 

the  expressions  f(n)  instead  of  f.  A  can  then  be  considered  as  a  modifi- 

! 

cation  of  the  MULT  algebra  (whenever  R  is  taken  to  be  the  SIMPLIFY  algebra) 

i 

9  !  ' 

having  as  its  generating  set  X  the  collection  of j maps  fn  such  that  f(Z-{Oj) 
fo}  and  as  its  basic  product  the  value  E  _  fjgj | instead  of  fngr .  However, 
if  A  is  considered  ate  an  algebra  over  the  polynomial  ring  defined  on  the 

f  1  I  1  r  y  ,  !  '  ’  ■  k 

set  (xj,  in  which  case  i_(a*  )fjr  =  afn_k  f°r  any  j  scalar  of  the  form  ax  anc 

i 

any  vector  f,  it  is  then  evident  that  the  algebra  is  a  modification  of  the 

j 

MULTILINEAR  algebra  with  generator  set  X,  whose  Scalar  ring  is  the  MULT 
algebra  with  generating  set  {x}.  In  either  case|  since  the  notation  being 

l 

used  to  represent  our  series  will  most  likely  differ  from  that  of  the  elemt 
f(n)  of  A,  a  short  input-output  package  will  have  to  be  supplied  to  conver' 
from  one  system  of  notation  to  the  ether.  But  this  is  always  the  case,  of 
course,  whenever  we  deal  with  isomorphic  algebras. 


Assuming  that  the  standard  differentiation  map  -j~-  :  fn  —  (n  +  l)/n  +  i 
is  associated  with  A,  it  then  follows  that  the  algebra  will  provide  a 
direct  procedure  for  reducing  ordinary  differential  equations  (in  x) 
to  their  difference  equation  components.  Moreover,  if  A  and  MULTILINEAR 
are  used  in  conjunction  with  one  another  for  handling  these  equations, 
first  transforming  them  and  then  reducing  them,  we  shall  find  that  we 
have  developed  a  highly  useful  procedure  for  examining  their  underlying 

structure.  .  . . „ _ _ _ _ , _ _ _ ^ _ _ 

Differential  Forms 

Let  R  be  the  MULTILINEAR  algebra,  having  as  its  generating  set  a 
sequence  of  the  form  {yi,...,ya},  the  choice  of  which  may  vary  with  the 
problem  under  consideration,  and  as  its  scalar  ring  the  MULT  algebra  with 
generator  set  {x1,...,xn].  Assuming  also  that  the  functions  being 
associated  with  MIJLT  and  MULTILINEAR  are  differentiable,  we  then  con¬ 
sider  the  exterior  algebra  A  of  the  free  module  M  generated  by  the  set 
[dxx  , . . .  ,dx„  }  over  R,  where  {dxj  , . . .  ,  dx,,  }  is  a  collection  of  n  atoms  or 
indexed  variables  (  ee  reference  2). 

Noting  that  the  non-scalar  elements  of  this  algebra  are  of  the  form 
uj  s  E«  v<1.  .  .  a,  ...i  dx,  -....dx,  (r  may  vary),  it  follows  that  the 

s  s  Ij,  *1  r  1  *  r  • 

algebra  is  a  model  of  the  algebra  of  differential  forms  for  some  n-manifold 

(see  reference  7).  This  in  turn  implies,  of  course,  that  the  differential 

n  da  > 

d(a)  ■  ^E  dxt  of  an  expression  a  in  R  is  an  element  of  this  algebra. 
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We  now  extend  the  definition  of  this  differentiation  mapping  in  the 
standard  manner  so  as  to  have  the  value  d(n)  =  £  d(at  ,.t  ) 

i  j  <  •  •  •  <1  r 

*dx. ,  *  . .  .-'dx.  for  any  vector  of  the  form  cu  =  E  ai,  ...i  dx, . 

1  'r  i1<...<ir  1  r  ll 

«...Adx,  ,  in  which  case: 

r 

(1)  d(otu1  +  3uj2)  =  odGuj)  +  3^(012)  for  any  vectors  and  uu2.  and 
any  elements  a  and  3  which  are  in  the  SIMPLIFY  ring. 

(2)  d(d(u>))  =  0  for  any  vector  uu  of  class  cr  (r  s  2). 

(3)  d((JU2A(u2)  =  d((jU2)Atu2  +  (“l)ruuiAd(ui2)  for  any  vector  aa  of  the  form 

iBj  =  E  at  .,(1  dxj  A...»dxlf  (r  fixed)  and  any  arbitrary  vector 

This  algebra  is  an  extremely  important  application  of  the  MULTILINEAR 
ring  R,  allowing  us  not  only  to  be  able  to  transform  many  highly  complex 
partial  differential  equations  in  a  systematic,  yet  simple  manner,  but 
to  transform  and  in  many  cases  reduce  multiple  integrals  also.  Moreover, 
if  the  collection  of  algebraic  homomorphisms  of  the  form  /: A  -•  A, where 
/( M)  .  M,  are  associated  with  A,  we  then  have  a  direct  procedur.  for  studying 
the  collection  of  matrices  associated  with  the  corresponding  restricted 
1  inear- transformations--/  |  M:M  —  M. - - — - — - - — — - - - 


On  the  Nature  of  Symbol  Manipulation  Languages 

During  the  last  few  years  a  number  of  studies  have  been  made  on  the 
structure  required  for  higher-level  languages  designed  primarily  for  an 
analyst  desiring  to  manipulate  various  mathematical  expressions  in  a 
variety  of  ways  (e.g.,  See  reference  3).  From  these  studies  it  has  become 
evident  that  the  general  list  processing  languages  (such  as  LISP)  are 
not  appropriate  for  processing  complex  mathematical  data.  This  is  due 
not  only  to  the  extreme  complexity  of  the  coding  that  necessarily  ensues 
in  the  analysis  of  an  intricate  problem,  but  also  to  the  fact  that  the 
input-output  procedures  most  ideally  suited  for  these  languages  are 
totally  inadequate  for  the  individual  primarily  interested  in  the  analysis 
of  a  problem  and  not  the  underlying  mathematical  and  input-output 
structure  involved.  On  the  other  hand,  it  has  become  equally  clear  that 
in  order  for  such  a  symbol  manipulation  language  to  be  constructed,  a 
list  processing  language,  or  equivalently  a  set  of  list  processing  routini 
will  necessarily  have  to  be  employed,  along  with  a  highly  sophisticated 
input-output  procedure  for  allowing  the  analyst  to  communicate  with  the 
machine  and  a  mathematical  structure  which  would  be  broad  enough  in  scope 
to  interpret  a  diverse  collection  of  expressions  in  a  fairly  simple  manne: 
We  shall  now  consider  briefly  the  underlying  mathematical  structure  that 
any  such  language  would  need. 
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In  general,  any  symbol  manipulation  language  must  always  contain 
at  least  two  modes  of  computation,  namely  the  ring  of  rational  numbers 
and  a  polynomial  algebra  over  a  set  of  variables  using  only  numerical 
coefficients  (in  which  case,  we  have  just  an  application  of  the  SIMPLIFY 
algebra,  of  course).  However,  because  of  the  diversity  of  the  mathematical 
concepts  encountered  and  the  techniques  that  are  needed,  these  modes  will 
obviously  never  suffice.  On  the  other  hand,  since  SIMPLIFY  was  found  to 
be  a  very  simple  application  of  an  extremely  general  model,  namely 
ALGEBRA,  and  such  a  model  in  some  form  will  always  be  needed  for  defining 
this  algebra,  it  is  evident  that  one  can  implement  the  model  itself 
so  that  we  then  have  at  our  disposal  not  only  the  SIMPLIFY  algebra,  but 
a  host  of  other  mathematical  procedures  also.  This  being  the  case,  an 
experimental  programming  language  called  FLAP  (for  FORTRAN-1 ike  algebraic 
list  processor)  was  developed,  using  the  ALGEBRA  system  as  its  underlying 
structure. 

The  FLAP  language  is  a  higher-level  language  (written  in  LISP)  which 
contains  not  only  the  model  ALGEBRA  but  also  a  complete  system  of  input- 
output  routines  for  intrepreting  the  data  under  consideration.  The 
language  contains  many  computational  modes,  including  the  standard  fixed 
and  floating  point  arithmetic  and  the  SIMPLIFY,  MULT,  and  MULTILINEAR 
algebras.  Its  format  is  similar  to  that  of  the  FORTRAN  system,  in  that 
it  makes  use  of  the  same  notation  of  the  statements  and  formulas,  the 
major  distinction  being  the  existence  of  a  mathematical  mode  setting 
statement,  appropriately  called  the  MODE  statement.  This  statement, 
which  has  for  its  argument  the  name  of  an  algebra,  in  effect  tells  the 
interpreter  that  until  another  mode  statement  is  made,  every  statement 
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is  to  be  regarded  as  being  dominated  by  this  algebra.  If  however,  no 

G 

MODE  statements  are  made,  the  standard  fixed  and  floating  point 
arithmetic  will  prevail  throughout  the  program  so  that,  in  effect,  we 
then  have  a  FORTRAN-1 ike  program. 

The  language  FLAP,  though  still  in  the  experimental  stage,  has 
already  been  found  to  be  extremely  useful.  This  is  due  not  only  to  the 
list-processing  capability  which  it  inherits  from  LISP,  but  also  to  the 
fact  that  it  has  been  formulated  on  a  general  analytical  structure 
(namely  ALGEBRAJ)  rather  than  a  set  of  distinct  computational  procedures. 

If 

Using  the  model  we  then  have  the  capability  of  being  able  to  easily  define 
any  other  computational  mode  (by  setting  the  arguments  X,  r,  v,  AXIOM, 
etc.}  as  well  as  being  able  to  modify  any  existing  ones. 


EXAMPLE: 


Suppose  that  we  wished  to  write  a  FLAP,  program  for  computing  the 

expression  —  [ay • sin(xy)] ,  where  y  is  a  variable  depending  on  the 
dxi  j 


variable  x  and  Ip.  is  a  constant.  Since  we  would  then  be  using  the 
MULTILINEAR  algebra  with  the  generator  set  {y}  and  the  scalar  ring 
MULT  (which  would,  in  turn,  have  the  generjator  set  {x}),  the  program 
could  be  written  as  follows:  r 
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CALCULATE 


INPUT  A,X,Y 
DUMMY  K 

MODE  MULTILINEAR  (X) , (Y) 

K  =  DIFF (A*Y*SIN (X*Y) , X) 

RETURN  K 
END 

CALCULATE (A, X,Y) 

The  following  value  would  then  be  printed: 

A*X*Y*COS (X*Y)*DIFF (Y, X)+A*Y**2*C0S (X*Y)+A*S IN (X*Y)*DIFF (Y,X) 

The  three  statement  cards  beginning  with  the  MODE  statement  contain  the 
instructions  for  computing  and  printing  the  expression;  the  other  cards 
are  used  only  for  declaring  and  binding  the  variables  concerned. 
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